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1. Introduction 

The connection between quantum dissipation and decoherence is a topic of longstanding 
interest p^EiEllHEiEiElEllS]. The main systems analyzed in this perspective are the 
damped harmonic oscillator, two level systems and quantum Brownian motion. For 
such systems the Hamiltonian description is not appropriate, and the most successful 
results come from the reduced description of a particle interacting with some type of 
reservoir [TO] . Classical understanding of the phenomenon is well-established, relying on 
Langevin or Fokker-Planck equations obtained by considering a particle interacting with 
a bath of independent oscillators [11]. The quantum counterpart of classical Brownian 
motion, however, has only recently been cast into standard equations [T2l H3] . 

Several approaches have been followed in order to obtain a quantum description of 
the dynamics of the Brownian particle: 

• a model-reservoir approach, leading to the famous Caldeira and Leggett master 
equation, which assumes the particle to be coupled to an environment described 
by a collection of simple harmonic oscillators [HE], and into which suitable terms 
can be added in order to produce a satisfactory Markovian equation of the required 
Lindblad form [31 H] , 
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• a dynamical approach modelling more closely the random collisions between the 
Brownian particle and the particles that make up its surrounding environment 

mum, 

• a measurement-based approach focusing on the information about the Brownian 
particle carried away by the particles of the medium during the collisions. This 
information is available, at least in principle, by monitoring the environment 
particles P El H21 H3] . 

The property of complete positivity, to be satisfied by a master equation for 
the reduced density operator of the particle, is a useful and stringent requirement in 
the study of subdynamics in quantum mechanics [14] . The various approaches are 
described in [12] . where the results obtained with the different methods are discussed. 
Of particular relevance is whether or not the proposed master equations are Markovian 
and of Lindblad form. 

Here the approach that we use to describe the quantum Brownian particle 
dynamics is the measurement-based one found in [T2], [131 US]. The collisions with the 
surrounding particles are considered to perform a random sequence of measurements 
feeding information about the position and momentum of the Brownian particle into 
the environment. Even if both of these quantities cannot be known with total precision 
at the same time, it is possible to simultaneously measure position and momentum 
by introducing some degree of imprecision for both. Using non-quantum-limited 
measurement techniques to represent the acquisition of this information [16] has led 
to the following master equation, in the limit of frequent collisions which make very 
weak joint measurements of position and momentum [12J: 



dp i 



p 1 * 

p 



[x, {P, p}] ~ [ £ > [ £ > Pi) - ^ [P> [P> P]] ■ (!; 



_2M' 

This equation is of the required Lindblad form provided that 

D PP D XX > (ft 7 /4) 2 . (2) 

Here M is the mass of the Brownian particle, 7 is the damping coefficient while D pp and 
D xx are diffusion coefficients given by 



D pp 7 

and 



Mk B T + ^(Kpf 



(3) 



8(A^) 



S- = ^T. (4) 

The particles forming the environment have mass m, R is the average rate of collisions 
and A a p and A a x represent the increase in the standard deviations due to the 
measurements of position and momentum over and above the intrinsic variances [T2] . 

Satisfaction of the condition of Eq. (|2J) using Eq. (jSJ) and (jlj) ensures that the master 
equation is of Lindblad form. The existence of analogous master equations obtained by 
other approaches, but with different expressions for D pp and D xx [HI HO], has been 
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described in [TS]. The Caldeira-Leggett master equation [TJ, is Markovian but is not 
of Lindblad form and this has been shown to lead to serious difficulties including, in 
particular, negative probabilities [TF1 [181 03] • These problems do not arise for Eq. (TjQ) 
because of the presence in D pp of two new terms other than the temperature dependent 
one, •yMksT. These depend on the variances A a p and A a x, and in particular the 
presence of a new double commutator term representing a position diffusion regulated 
by D xx . The origin of these terms is the inherent spreading in position(momentum) 
which occurs when a measurement of momentum(position) is made. This point has 
been discussed further in [13]. 

The main aim of this paper is an analysis of the decoherence dynamics associated 
with Eq. ([Tj) . In particular we find the exact solution of the master equation, and then 
use this to illuminate the role of both the D pp and the hitherto largely unconsidered D xx 
terms in dissipation and decoherence. To this end we will consider two initial states for 
the Brownian particle: a single Gaussian wave packet and Schrodinger cat state. By way 
of comparison we also consider cases where D pp and D xx , can be varied independently, 
which can, for example, furnish the Caldeira-Leggett dynamics when D pp = jMksT 
and D xx = 0. 

The paper is organized as follows. In Sec. [2]we solve the master equation, providing 
a general scheme for obtaining time evolution from general initial conditions. In Sec. [3] 
we apply this scheme to two initial configurations, and discuss their dynamical evolution. 
In Sec. H]we consider briefly what happens when the Lindblad condition given by Eq. (j2J) 
is not satisfied. In Sec. [5] we summarize and discuss our results. In |Appendix A| and 



Appendix B we collect some of the lengthier calculations. 



2. Solution of the master equation 



In this section we solve the master equation by introducing a characteristic function. 
This procedure has been previously used to obtain a formal solution of the Caldeira- 
Leggett master equation, and is described, for example, in [19]. The main difference 
here is that there is an additional term depending on the position diffusion D xx . 
In the position representation Eq. ([1]) takes the form: 



dp(x, x', t) 
dt 



" ih f 


Q2 


d 2 \ 


2M \ 


dx 2 


dx' 2 ) 


D xx 


( d 2 
\ dx 2 


d 2 
dx' 2 




D 



PP I t\2 
' X — X ) 



— p(x,x',t) 



(5) 



where p(x,x',t) = (x\p\x'). This second order linear partial differential equation can be 
greatly simplified by moving to a (k, A t ) representation [20] based on introducing the 
characteristic function p(k,A t ,t): 

hA t " 



p(k, A t , t) = tr(Dp) = tr(exp[z(/cx + A t p)]p) 



dx e 



p 



' hA t 

x H , x 



t . (6) 
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In this new representation we obtain a first order partial differential equation in the 
form 



dp(k,A t ,t) 
dt 



k d d 
-j A t - 



DppA t D xx k 



p(k, A t ,t) 



(7) 



MdA t 1 "dA t 

In Appendix A the method of characteristics [21] is used to solve this partial differential 
equation exactly. The solution is 

k 



P (k, A t ,t) = p[ k, A t (i-r 

D 



jM 



exp 



where Ar 



D xx 



'pp 



t + 



D 



pp 



M 2 7 3 



r,o x 



^ pp ^ 2 -kA t — ^PpTj'Z £) ^2 



Mr 



2l 



(8) 



A t (l - T) + -^r and T = 1 - exp(- 7 t). 



^0 — "tv- - i i lM > 

We can use Eq. ([6]) and its inverse to move to and from the coordinate and 
(k, A t ) representations at will, and Eq. (jS|) to obtain the time evolution in the (k, A t ) 
representation. This allows us to compute the time evolution of a general initial state 
in the position representation by following the procedure: 

p{x, x', 0) -> p(k, A , 0) -> p(k, A t , t) -> p(x, x', t) . (9) 

In the next section we will use this procedure to solve the master equation for 
two initial states which can be written in the coordinate representation as a sum of 
exponential terms of the form 

p(x, x',0)— exp — A (x — x') 2 — iB (x — x')(x + x') 

-C {x + x') 2 - iD (x - x') - E (x + x') - F ] . (10) 

We here apply the scheme outlined by Eq. ([9]) to obtain the time evolution for states 
of this kind. Using Eqs. (EJ) and (El) we find the characteristic function at time t in the 
(k, A t ) representation to be 



p(k, A t , t) = exp -a t k 2 - ib t kA t - c t A 2 t - id t k - ie t A t - f t 



where the various coefficients follow the time evolution given by 

2 / 



a t = a + b —— + c 

M7 



Mi. 



+ \D XX + 



D 



pp 



, / s 2r(i-r) d vv t 2 

b t =b (l-T)+c —^7 L + " 



M 2 7 2 
c t = 



t 



D 



pp 



M 2 -f 3 



d t = d + e 



MY 



M7 ' M7 2 
e i = e (l-r), f t = fo 



co(i-r) 



D PP T(2 - r) 



27 



(12) 



and where the relation between small and capital coefficients is given by 

1 . B AAqCq + B 2 2 Eq 

n, c = — n 



a 
e 



16C ' 4C 
2C -Do — £q-Eo 



4C 



rf 



2C n 



% exp(-Jo) = exp(-Fo) exp 



4C 
E 2 N 



7T 



4Cn / 2 V C, 



7T- (13) 
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Transformation back to the coordinate representation Eq. fill I) provides the density 



For density matrices which correspond to superpositions of states centred on 
different positions, each of which is individually of the form given by Eq. ( 11 01) . the 
linearity of Eq. (J7J) means that it is still possible to follow the procedure described by 
Eq. (J9]) and so to compute the time evolution. 

3. Applications 

In this section we consider two physically interesting initial states, the simple Gaussian 
wavepacket and a superposition of two such states, which forms a Schrodinger cat state. 
Using the procedure described in the preceding section to compute their time evolution, 
we focus on the role of the momentum and position diffusion terms, proportional to D pp 
and D xx respectively. In particular many of our results will be expressed in terms of 
these coefficients and will not depend on their explicit form [Eqs. (J3]) and (Jl])] in terms of 
the physical parameters of the system. These results have a general validity, therefore 
for any master equation with the same form as Eq. (JT]) irrespective of the sizes of the 
diffusion terms. Thus our method does not just provide the solution to Eq. (jl|), but also 
an infinity of other master equations whose diffusion terms do not necessarily satisfy 
the Lindblad condition (Eq. (J2])). One particular example is the Caldeira-Leggett master 
equation obtained by putting D pp = ^MhsT and D xx = 0. 

3.1. Single Gaussian wave packet 

Consider a minimum uncertainty Gaussian wave packet centered at position x = 
and momentum po, with initial spreads in position and momentum Axq and Ap which 
satisfy the uncertainty principle Ax Ap = |, 



matrix 



p(x,x',t)= exp — A t (x — x') 2 — iB t (x — x')(x + x') 

-C t (x + x'f - iD t (x - x') - E t (x + x') - F t 
with the following relationships 



(14) 




(15) 



(16) 



p(x,x',0) 




1 



exp 




+ i 



.Po(x — x') 

1 % 



(17) 
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By comparing this initial reduced density matrix with Eq. (HOp at t = we can identify 
the required coefficients as 

^o = C = -^, B = E = 0, A) = -x> e M-Fo) = ^==- (18) 

By using Eq. ([6]) we can obtain the corresponding initial condition in the (k, A t ) 
representation, which has the exponential form of Eq. (Hip for t = with 

«o = — b = 0, c = — d = 0, e = -p , /o = . (19) 



^From Eq. (fl2l) we next compute the time evolution of the various coefficients, and 
then return to the coordinate representation using Eqs. (|14p and (|15p . Thus we obtain 
the spatial reduced density matrix at time t, from which it is also possible to obtain 
the corresponding density matrix in the momentum representation by double Fourier 
transformation. The solution provides a simple means of obtaining the time evolution 
of the average of x and p, and of their variances shown below and plotted in Fig. 1 for 
physically reasonable parameter values, 

(x) t = - dt, (P)t = ~et, Ax 2 = 2a t , Ap 2 t = 2c t . (20) 

These equations could equally well be rewritten in terms of the capitalized coefficients 
using Eq. (1161) . but the evolution is given most simply in terms of the initial conditions 
using Eq. (}T2l) . The expectation values become 



= 77 ' = Poexp(-7t), Ax t = Ax + 



1 — exp(— 7t) 



7 



2 

+ 



-pp -PP 



3 



'1 - exp(- 7 t)) 2 1 



2 



Ap? = Ap* exp(-2 7 t) + D pp [1 6XP( 27t)] . (21) 

7 

As would be expected, the diffusions D pp and -D^. do not affect the mean values of 
the position and momentum. Of particular interest, however, is the evolution of the 
variances, whose dependence on D pp and D xx , can be simply found from Eq. (1211) for 
both short times (t <C 7 _1 ) 

Ax 2 « Ax 2 + 2D xx t, Ap\ « Apg + 2 (l^ - 7 Apg) t , (22) 



and long times (t 7 



A ^ w2 ( I, - + ]^) t ' ^"^p (23) 

Eqs. (122 p and (I23p reveal the critical role of D pp and -D^. In particular, for small times 
Ax 2 depends only on D xx , which shows that the D xx term has an important role in this 
time region. 

It is possible to calculate the evolution of the purity, ti(p 2 ) t , of the initial state of 
Eq. ([IT]) . This quantity, bounded by and 1, is related to the linear entropy and is 
equal to 1 for pure states. Any difference from 1 means a loss of purity of the state. The 
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<x> t <p>, Ax t Ap t 




Figure 1. Figures on the left: evolution of the average of x and p. Figures on 
the right: evolution of the variances Axt and Ap t . Parameter values (in SI units): 
M = 5.01 x 10~ 22 , m= 5.01 xl0~ 26 , k B = 1.38 x 10~ 23 , h = 1.06 x 10~ 34 , 7 = 
1000, R = jM/2m = 5 x 10 6 , T = 300K, p = 5.01 x 10~ 26 , Ax Q = 
0.73 x 10~ 7 , A Po = 7.26 x 10~ 28 (A CT p) 2 = (A^x) 2 = | X n, u = 10 4 . These 
parameters remain the same in all the following figures. 

tr(p 2 ) t 

10- 17 \ 




Figure 2. Behavior in time of tr(p 2 )t. 



purity of the Gaussian state is plotted in Fig. 2 as a function of time. The figure shows 
that the decoherence process occurs on a time scale much smaller than the relaxation 
time of the particle, which is of the order of 7 -1 . The dynamics of our system are 
then described by the particle density matrix time evolution as a rapid transformation 
from the pure initial state (TT7T) into a statistical mixture, as is shown in Fig. 3 in the 
coordinate representation. 

To investigate further the dynamics of this loss of coherence both in x and p space 
we can compute the spreads M. 2 x (t) and Mp(t), and coherence lengths, £ x (t) and £ p (t), 
using their general definition in [22] in terms of traces over the density matrix of the 
particle 

Tr(p 2 x 2 ) + Ti(pxpx) 




M 2 x (t) = ^ * - 2 ( 4^ ) , (24) 



m - ^t^im , (25) 

with similar expressions for the momentum spread and coherence length. For pure states 
p 2 = p and therefore both spread and coherence length become equal to the respective 
width of the state, e.g. M. x {t) = C x (t) = Ax t . However, in presence of the interaction 
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Figure 3. Left: Absolute value of the initial density matrix (|17p . Right: \p(x, x' , t*)\, 
where t* <C 7 _1 . The suppression of the off-diagonal terms reflects the onset of 
decoherence before significant spreading has occurred. 



the state of the particle loses its initial purity and the two quantities differ. While M. x {t) 
gives the extension of the state, C x {t) gives the zone inside the state extension, where 
coherence has not yet been lost at time t |19]. 

For an initial Gaussian wave packet the spread A4 x (t) corresponds to the width of 
the reduced density matrix along the main diagonal, M. x {t) = Ax t = J p(x,x,t)x 2 — 
(/ p(x, x, t)x) 2 , while the coherence length C x (t) gives analogously the width of the 
reduced density matrix along the main skew diagonal C x (t) = J p(x, — x,t)x 2 — 
(fp(x, -x,t)x) 2 . 

The ratio C x (t) / A4 x (t) gives a dimensionless measurement of the loss of coherence. 
It is interesting to note that for an initial Gaussian wave packet [23], this ratio is directly 
connected to tr(p 2 ) and is equal in both the position and momentum representations. 
This property is found also in our system where the following relations are obtained 

tr(A = ! _ 5 „„ = M> = M> = V 2 . ( 26) 
Ax, Ap, 

^From Eqs. (12^|) and (|2TJ|) we obtain 



£lW A P « = £ t(t )A Il = ^_L_>| (27) 

which can be seen as a particular case of the generalized uncertainty relation derived in 
[22]. 

The squares of the coherence lengths in the two representations are given by 

^iS^? md *'^iS- (28) 

It is instructive to consider the evolution of the coherence lengths at short times 
(t « 7" 1 ) 

D' 



C 2 x (t) « Ax 2 



Dpp 7 ) t 



Cl(t)*A P Ul-2^t , (29) 
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and for large times (t 7 ^ 

ft 2 



h 2 1 
4D. 



4Ax 2 



(30) 



Eqs. ( l29l and ( !30l show the role of D pp and D xx in these evolutions, i.e. in the 
decoherence process. In particular, for small times the momentum coherence length 
Cp(t) depends only on the position diffusion D xx , showing again the relevance of this 
term in this time region. This behaviour of C 2 (t) can be understood by comparing it 
with that of £ 2 (t) which, as it is independent of D xx , holds also in the Caldeira-Leggett 
model. Indeed, from Eqs. (1291) and ( 1221) we see that, as C 2 x {t) and Apf both depend on 
the factor (D pp / Ap% — 7), and the same occurs for C 2 {t) and Ax 2 , which both depend 
on Drr. 



3.2. Schrddinger cat state 

The second configuration that we investigate is an initial Schrddinger cat state with a 
model wavefunction of the form 

1 



1>(x, 0) 



x 



exp 



x 



1 + exp 

2 



P 



4a 2 



Mv 
— i——x 
n 



(27nr 2 )-i 
+exp 



x + 2 



4a 2 



Mv 
+ i——x 
n 



(31) 



where a is the width of two wave packets initially placed at a distance I with one 
moving towards the other with an initial velocity v. Such states are interesting because 
of their potentially long-range coherence properties and the extreme sensitivity of this 
to environmental decoherence [24] . 

In the absence of any interactions the behaviour of the diagonal reduced density 
matrix elements is given by 



p(x, x, t) 



1 



2 1 



2tt a 2 + 



4M 2 a 2 



exp 



X + \ - vtf 



2 a 2 + 



h 2 t 2 

4M 2 a 2 



+ exp 



X 2 


+ G 


\ 2 
-vtj 


2( 


a 2 + 


h 2 t 2 \ 


4M 2 a 2 ) 



exp 



cos 



[X 



\ + vtf 



2 a 2 + 



+ 



2.2 



i, 2 t 



AMva 2 

h 



4M 2 a 2 



hit 

2a 2 M 



+ 



2 a 2 + 



h 2 t 2 

4M 2 a 2 



X 



(32) 



This probability distribution is the sum of three contributions. The first two clearly 
correspond to a pair of separately expanding (freely spreading) wave packets, while 
the third term represents an interference term, responsible for the central peak present 
in Fig. 4, which is a plot of Eq. ( 1321) . This is exactly as would be expected for such 
undamped evolution. 



In |Appendix B| the time evolution is computed from the initial state of Eq. (1311) 
using both the procedure described in Eq. (jUJ) and the linearity of Eq. (J7J). Along the 
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x = x' diagonal we find 

i 

p(x,x,t) = 



1 + exp( 



+ exp 



a:' 



- \d 



T |2 

t 



4a, 



— i ex P 

_ 2MW 
a 2 



Ad t 
I 2 



+ exp 



(x - d t f 
Ad t 



8a 2 



cos • 



\dt\x 
~2aZ 



(33) 



with the evolution of the various coefficients given in Eq. (IB. 4ft . Eq. (1331) reduces to 
Eq. (1321) in the absence of any interaction, as may readily be verified by substituting 

: 0. This probability distribution is again the sum of 



t and D pp 



lim 7 _ r/7 

three contributions. In order to investigate the behavior of the interference term, we 
compute the attenuation coefficient W t [25], defined as the ratio of the factor multiplying 
the cosine interference term to twice the geometric mean of the first two terms: 



W t = exp 



d\ + \d t 



2M 2 v 2 a 2 



(34) 



4a t 8a 2 h z 

In the free case Wt = 1 for all times, corresponding to a full coherence. In the presence 
of interaction W t decays quickly, as pictured in Fig. 5, which shows the rapid destruction 
of the interference with time. 

By expanding W t , for small times, in a Taylor series we find: 



4cr 



+ 



AM 2 v 2 a r 



Tr 



t. 



(35) 



This equation shows that for very short times the attenuation factor decays with a 
characteristic time td given by 



r 2 



AM 2 v 2 a 2 ^ 



m= ^l^ + ^-F-Jj • (36) 

The interaction with the environment leads to the destruction of the interference term. 
The initial decay of Wt, which characterizes this, is due solely to the presence of the 
D xx term in Eq. (TO . All of the decoherence occurs on a very short timescale and is due 
entirely to this diffusion. 



Frictional quantum decoherence 



11 





Figure 5. Evolution of the attenuation coefficient Wt- 



In order to investigate further the dynamics of this loss of coherence in x and p 
space we compute the spreads, M^{t) and Ai^it), and coherence lengths, C x (t) and 
C p (t) defined in Eqs. ((21]) and (1251) . By using Eq. (IB .71) we obtain the spatial spread and 
the coherence length at time t 

1 
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2F 4 + exp 



cos a 

'A 2 
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2F t + exp 
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where a 



[C t D t D t + A t E t E t ) /AA t C t + F t + F t and /? = [C t D* - (7* A + ^A - 
/8AtCt- In order to obtain the same quantities in momentum space we need the 
corresponding reduced density matrix. If we perform a double Fourier transform we 
obtain p(p,p',t) from p(x,x',t). The result has the same form as Eq. ( IB. 71) with the 
following substitutions: 

B t D t + 2A t E t 



A^ 



E t -*2h 



B t E t 



-B t 
Z t 
2C t D t 



Z t 



a 



E t 



Ct D 



>2h- 



2h 



B t E t - 2C t D t 
Z t 



z t 

Ff—*Ft- 



, A 

lnZ t 



+ 
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2 CtD 2 t -B t D t Et-A t E 2 t 



lnZ f 

F t -^F t +— -+Ah 2 



C t D 2 t - B t b t Et 



-2K- 



AtE} 



B t D t + 2A t E t 



(39) 
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where Z t = ^{B 2 + 4A t C t )h 2 . By using these substitutions in Eqs. (157j) and (1551) one 
obtains the corresponding -Mp(t) and 

These quantities satisfy .M^O) = C x (0) and -M p (0) = C p (0) at the initial time, 
while for large times their increase is the same as the corresponding quantities Ax t , 
l x (t), Ap t and l p {t) in Eqs. ( I2"5j) and (15U1) found for the single Gaussian wave packet. In 
particular it has been shown that a generalized uncertainty relation holds [22] 

C x (t)M p {t) > 1. (40) 

This product is plotted in Fig. 6 for our cat state, which shows how the uncertainty 
relation is satisfied at all times. 

Xx(t)M P (t) 




2 



Figure 6. The evolution of the product £ x (t)A4 p (t). 



4. Unphysical parameter region 

As was stated in the introduction, one of the key differences between Eq. ([T|) and the 
Caldeira Leggett equation is the presence of the D xx term. The latter master equation, 
of course, is not of the Lindblad type and pathological behavior has been observed by 
several authors [T9J, [13J, [18] if 7 is too large or the wave packet width too small (smaller 
than the thermal de Broglie wavelength of the object). If we consider the time derivative 
of the linear entropy for small times, then we would expect a positive value for an initial 
pure state, for which Su n = 0. In our model we find this time derivative to be 

^-Si in \ t=0 = --^-trp 2 | t=0 = -2(p) = -^Ap 2 D xx + ^Ax 2 D pp - 7 , (41) 
at at n n 

where the brackets indicate the average over a general state of the Brownian particle 
and we have used Eq. (pQ) to obtain the last equality. If we use AxoAp = |, then the 
positivity of Eq. (f4~Tl) is assured for all possible states only if 

D XX D PP > fa) . (42) 



This is the same condition as that found in [13] by requiring an initial reduction in the 
probability of remaining in the initial pure state. 
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Even if in our model this condition is satisfied by our parameters of Eqs. and 
(JH) it is interesting to work near the region of its validity. In fact we can use our master 
equation Eq. (0Q), and look for anomalous features when Eq. ( )42l) is not satisfied. For 
example, after substituting D pp = ; we can var Y Q around one. In Fig. 7 Eq. (HTI) 

is plotted as a function of q. 




Figure 7. Behavior of the derivative of the linear entropy at t = as a function of q. 



This figure clearly shows that the rate of change of linear entropy is negative at 
small times if the Lindblad condition is not satisfied. 

5. Conclusions and discussion 

In this paper we have used the Lindblad master equation for Brownian motion found 
in [121 13] to analyze wave packet dynamics. To this end we have provided a simple 
and clear scheme by which we can obtain the exact time evolution starting from general 
initial conditions. 

We have applied this procedure to find the time evolution of two initial states: the 
first in which the Brownian particle is represented by a Gaussian wave packet and the 
second, in which it is represented by a Schrodinger cat state. In each case we have 
provided expressions for the relevant quantities in both the position and momentum 
representations. There are complementary aspects of the two representations, such as 
generalized uncertainty relations not found when one focuses only on the spatial features. 

We have analyzed further the dissipation and decoherence, in particular focusing on 
the reduced density matrix evolution in the coordinate representation. We have obtained 
the evolution of the wavepacket widths and the coherence lengths for both initial 
conditions and in both the position and momentum representations. The Gaussian wave 
packet shows a very rapid decoherence on a timescale much shorter than the relaxation 
time of the system. In the momentum representation this decoherence depends only on 
the position diffusion coefficient, D xx , present in the master equation. The D xx term is 
also the only term responsible for the increase in the spatial width of the wave packet 
for small times. The term containing this coefficient is absent in many previous master 
equations which have been used to describe friction and Brownian motion. 
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The evolution of a Schrodinger cat state without any frictional effects shows 
coherent oscillations caused by interference between the two components of the wave 
function. Our analysis shows that such interference is damped on timescales again much 
shorter than the relaxation time; decoherence is very rapid. We have quantified the 
decoherence using an attenuation coefficient for the oscillatory terms. This attenuation 
coefficient also depends on D xx . The results of the computation of generalized variances 
and coherence lengths show that for large times these quantities behave in a similar way 
to the single Gaussian wave packet case, and that a generalized uncertainty relation 
between variances and coherence lengths holds for all times. 

Finally we have generalized the system to look at the case where the product of 
the diffusion coefficients is not large enough to guarantee that the master equation is 
of Lindblad form. Here we see a clear signature of unphysical behaviour in the linear 
entropy, which is associated with probabilities outside the physically-meaningful range 
between and 1. 

The measurement-based quantum description of friction illustrated in this paper 
provides a general framework for investigating the role that the various terms in the 
master equation play in decoherence. It is clear from our analysis that the two extra 
diffusion terms not associated with the temperature of the system are necessary to ensure 
complete positivity of the density operator at all times. This is consistent with previous 
work [12], [13]. The minimum sizes of these terms are governed by an uncertainty relation, 
in line with their wholly quantum origin. Diffusion in one observable is associated 
with localization in the complementary one. These localizations occur each time a 
measurement is made. The consequent diffusion, and in particular that of position, 
which has no analogue in the Caldeira-Leggett reservoir-based approach, must be taken 
account of in any complete quantum description of friction. The cost of not doing 
so is illustrated emphatically here; the resultant incomplete equation cannot describe 
decoherence correctly, because it is not valid on the short timescales during which 
decoherence occurs. 

The generality of the measurement-based approach reflects the generality of the 
Kraus formalism of quantum measurements on which it is based, which makes no 
reference to any particular measurement device. Consequently the theory presented 
here is not specific to any particular frictional or Brownian system. Such a linkage 
could in principle be found for particular systems, and would amount to an ab initio 
quantum theory of friction. No such theory is known at this time. 
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In this appendix we use the method of characteristics [2T] to solve the master equation 
of Eq. ((Tj) . The first step is to rewrite Eq. (J7J) as 

k \ dp(k,A t ,t) 



, dp(k,A t ,t) , ( 

1 m + [ jAt 



M 



DppAf + D xx k ) p(k, A t ,t). 



The curves in the t, A t plane, parameterized by I and defined by the relation 
dt dA t 
T " 



A. 



A. 

M 



dl 



(a.i; 



(A.2) 



are called the characteristic curves of the partial differential equations. Eq. ( 1A.2I) may 
be written as: 

dt dA t dp(k,A t ,t) 

T " " 



7A t - 1 



(A.3) 



(D pp A 2 t +D xx k 2 )p(k,A t ,t) ' 

This pair of equations, valid on each characteristic curve, enables a general solution 
of the partial differential equation flA.ll) to be found as follows. We perform a first 
integration using the first equality of Eq. ( 1A.3|) . finding for one arbitrary constant W: 

k \ 



W 



-it 



At 



from which it follows at t — 

k 



W = An 



7M 



(A.4) 



(A.5) 



as k is independent of time. Now we perform a second integration using the second 
equality in Eq. (1A.3I) . finding a second arbitrary constant Z: 

A t 



p(k,A t ,t) 



cxp 



D PP A 2 t 



D pp kAt 



27 



M7 2 



7M 



from which it follows at t — 



1 



p(fc,A ,0) 



exp 



D pp Aq D pp kA 



27 



M7 2 



Ac 



(A.6) 



(A.7) 



By using Eq. (jA~7Fj) in Eq. flA.6|> we find for p(A;, A t ,t): 



p(k, A t , t) = p(k, A , 0) exp 
A t 



D pp (A 2 



7M 



27 



D pp k (A t - A c 
M7 2 



a --4t 



7M 

In order to express A t as a function of A we use Eqs. ()A.4j) and flA.5D . finding 

fc 



A r 



A t e"^ + --fl- 
* 7M V 



-7* 



(Ai 



(A.9) 
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Substituting the previous expression for A in Eq. flA.80 . we finally obtain for p(k, A t ,t): 



P (k, A t , t) = P (k, A t (i - r) + -^r, o j x 



(A.10) 



exp 



D + Dpp 

^ XX ~ 



t + 



D PP (l 2 



M 2 7 2 J M 2 7 3 V 2 
where T = 1 — exp(— ft). 



M-f 2 



k 2 _ t^ kAf _ ^P r ( 2 - r ) A 2 



27 



t ( ' 



Appendix B. 



In this appendix we obtain the time evolution of the Schrodinger cat initial state given 
in Eq. ([31]) . 

Computing the reduced density matrix p(x, x',0) corresponding to Eq. (|3T|) and 
using Eq. §6§ to move to the (A;, A t ) representation we obtain: 

- 2 h 2 

1.2 " \ 2 ; ' ; , r .. \^ ( . 



p(fc,A ,0) 



1 



2{l + exp[-£ s ] 



exp 



-Y k2 --^ A ° + ll 2 k ~ tMvA{ 



+ 



exp 



a 



—k l 



h 2 



—A 2 - i-k + MvA 



+ exp 



<7 2 h 2 a9 2Mva 2 r hi 

— k ^An - k ^A 

2 8(x 2 h Aa 2 



+ exp 



a 2 , h 2 . o 2Mwcr 



8a 2 



2 , U a ' 
+ — Ao 



2 " 8a 2_u ' H " ' Aa 2 
The form of this initial condition is of the kind: 



x exp 



2MVa 2 I 2 



h 2 



8a 2 



(B.l) 



p{k,A ,0) = Y / pAk,A ,0), (B.2) 
where comparing with Eq. ( 11 II) we have for the various coefficients the initial values 



<7 



1,2,3,4 
«0 = yj "0 



1,2,3,4 



0, c 



1,2,3,4 



8a 



2 5 u — u 



1 d 3-^-^ 
, u — u — J 



e = ~ e = Mv ' e n = -en 



c3,4 



/, o — Cq 



^1 — exp 


" d 2 ' 






) 




8a 2 _ 





+ 



2MVa 2 I 2 



1 — exp 



( / 2 



8a 2 



Y 2 ' 8a 2 ' (R3) 
By using Eq. (fl~2l it is possible to compute the time evolution of all the coefficients of 
the four parts of Eq. (IB. 21) : 



1,2,3,4 

a t ' = a t = — + 



a 2 h 2 7 2 



2 8a 2 M 2 7 2 



M 2 7 2 



M 2 7 3 V 2 



7,1,2,3,4 7 

6* = h 



h 2 2r(i-r) ^ z^r 2 1AM 



^CT 2 M 7 



+ 



h 2 



Mr 



■(i-r) 2 + 



2 , D pp r(2 - r) 



dl = -d 2 = d t = - l - + v -^, dl = -dt = d t = -i. h > Aa2Mi) 



8a 2 y ~ 2 7 
(2Mva 2 hi T \ 



ft/ 



ej = -e 2 = 6* = - T), e? = -e* = e~, = -<— (1 - T), 



^1,2 j?l,2 (■ ^d,4 f 

Jt ~~ J ~~ /*> //I ~ JO ~ Jt ■ 



3,4 _ z-3,4 



(B.4) 
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Then using Eq. (1151) we can move to coordinate representation obtaining for the reduced 
density matrix 



4 4 

p(x, x', t) — Pj(x, x', t) = exp —A J t (x — x') 2 — iB\[x — x')(x + x') — 

Cl{x + x') 2 - zl^'(z - x') - £|(x + x') - F/] , 
where 

a 1,2,3,4 7 4a 4 Q — b 2 D l, 2,3,4 5 &i /^fl,2,3,4 ^ 1 

A = A t = — — 2Z , -t>t = tS t = --tt— , Ot = Of = — — , 

Ah a t 4:ha t 16a t 

, - 2dt6f — bfdt n 2dt6t — b+dt Q 

D 1 = -D, 2 = D t = — — ~ — — , D, 3 = -D4t = — — ~r — — = A, -Ey = -£? = Ft 

Ef = -Et = & = exp(-F>< 2 ) = exp(-F,) = ^^exp[-^" 

<Aa t zyVat [ 4a 4 

4cif 



exp(-F t ) = exp(-Ft) = - exp 

Next from the last equati 
p(x, x', t) = exp 
{exp 
+ exp 



ion in Eq. (IB. 51) it follows that 
— A t (x — x') 2 — iB t (x — x')(x + x') — C t (x + x') 
-iD t {x- x') - T ' ' ' " ' ' ' 



X 



exp 



+iD t (x - x') + E t (x + x') - F t 



(B.5) 



d t 

4a*' 



(B.6) 



x') -E t (x + x') ~F t 
-W t (x - x') -E t (x + x') - F t ] + exp[+W t {x - x') + E t (x + x') - F t ] } . (B.7) 
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